System size and centrality dependence of charged hadron transverse momentum 
spectra in Au+Au and Cu+Cu collisions at y/s NN = 62.4 and 200 GeV 
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We present transverse momentum distributions of charged hadrons produced in Cu+Cu collisions 
at y^sT/Tv" = 62.4 and 200 GeV. The spectra are measured for transverse momenta of 0.25 < pr < 
5.0 GeV/c at v /s7^~ = 62.4 GeV and 0.25 < p T < 7.0 GeV/c at ^/s^T/ = 200 GeV, in a pseudo- 
rapidity range of 0.2 < r\ < 1.4. The nuclear modification factor Raa is calculated relative to p+p 
data at both collision energies as a function of collision centrality. At a given collision energy and 
fractional cross-section, Raa is observed to be systematically larger in Cu+Cu collisions compared 
to Au+Au. However, for the same number of participating nucleons, Raa is essentially the same in 
both systems over the measured range of pr, in spite of the significantly different geometries of the 
Cu+Cu and Au+Au systems. 
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The yield of charged hadrons produced at mid-rapidity 
in Cu+Cu collisions at energies of ^/s^Tv" = 62.4 and 
200 GeV has been measured with the PHOBOS de- 
tector at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory. The data are 
presented as a function of transverse momentum (pt) 
and collision centrality. The goal of these measurements 
is to study the possible formation of a new form of matter 
through modification of particle production compared to 
nucleon-nucleon collisions at the same energy. 

This measurement was motivated by results from 
Au+Au collisions for ^/s7^ = 62.4, 130 and 200 GeV. 
Hadron production at these energies was found to be 
strongly suppressed relative to expectations based on an 
independent superposition of nucleon-nucleon collisions 
at p T of 2-10 GeV/c [1 II S @ . The modification of 
high-pT hadron yields has commonly been investigated 
using the nuclear modification factor, Raa, defined as 



Raa(pt) 



d 2 N A A/dp T dri 



(1) 



{Ncoii) d 2 a pp /dpTdr] 
A value of Raa = 1 is obtained if particle production 



scales with the average number of binary nucleon-nucleon 
collisions, (N co u), within a heavy-ion collision. Instead, 
for the production of charged hadrons in central Au+Au 
collisions at ^fs^ = 200 GeV, values of Raa ~ 0.2 are 
observed at p T = 4 GeV/c @|10. 

Such a suppression had been predicted to occur as 
a consequence of the energy loss of high-pr partons in 
the dense medium formed in Au+Au collisions 6j. This 
hypothesis is also supported by the observed absence of 
this effect in deuteron-gold collisions at the same collision 
energy 0, H H . 

The results presented here for Cu+Cu collisions at 
y/s NN = 62.4 and 200 GeV bridge the gap between 
the Au+Au and d+Au systems, allowing a unique 
examination of the dependence of high-py suppression on 
system size. A careful comparison of Cu+Cu and Au+Au 
spectra to model calculations at high-p^ can elucidate 
the dependence of absorption on path length, especially 
as one expects a different distribution of path lengths for 
jets produced in the two systems even for the same N part . 

The data were collected using the PHOBOS two- 
arm magnetic spectrometer. Details of the experimental 
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TABLE I: Details of the centrality classes used in this 
analysis. Bins are expressed in terms of percentage of the 
total inelastic Cu+Cu cross-section. 
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FIG. 1: Mean number of collisions per participant pair 
for Au+Au (open symbols) and Cu+Cu (filled symbols) at 
^s^T = 200 GeV. 

setup can be found in The primary event trigger 

used the time difference between signals in two sets of 10 
Cerenkov counters, located at 4.4 < \rj\ < 4.9, to select 
collisions that were close to the nominal vertex position 
along the beam-axis. 

For the analysis presented here, events were divided 
into centrality classes based on the total energy deposited 
in the octagon silicon detector, covering pseudo-rapidities 
< 3.2. A full detector simulation using HIJING 
events 0, was used to estimate (N part ) for each 
centrality class, and the corresponding (N co u) values 
were obtained from a Monte Carlo Glauber calculation 
[T2I FLU ]. For these calculations, as well as for the 
determination of Raa at 62.4 and 200 GeV. we used 
ap n p el = 36 ± 1 and 42 ± 1 mb, respectively 15]. The 
results are listed in Table I. The systematic uncertainty 
on these values comes primarily from the uncertainty 
on our estimate of the measured fraction of the total 
inelastic cross-section, which is calculated by a variety of 
methods in a range of pseudo-rapidity regions (lif. The 
uncertainty also covers the ran ge o f efficiencies measured 
using HIJING and AMPT 0113. 

Fig- IH demonstrates that the value of (N co u) is essen- 
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FIG. 2: Top: Invariant yields for charged hadrons from 
Cu+Cu collisions at ^/s NN = 200 GeV, in the pseudo-rapidity 
interval 0.2<?7< 1.4 as a function of pr for 6 centrality 
bins. Bottom: The same for ^/s NN — 62.4 GeV. For clarity, 
consecutive bins are scaled by factors of 5. Statistical and 
systematic uncertainties are smaller than the symbol size. 



tially the same in Au+Au and Cu+Cu for the same value 
of (Np ar t). Thus, the comparison of these two systems 
does not allow one to distinguish between scaling with 
participants or collisions. 

The event selection and track reconstruction procedure 
for this analysis closely followed the previously published 
analysis of Au+Au collisions at y/s NN — 62.4 GeV y|. 
Events with a primary vertex position within ±10 cm of 
the nominal vertex position were selected. Only particles 
traversing a full spectrometer arm were included in the 
analysis, resulting in a low transverse momentum cutoff 
at p T w 0.2 GeV/c. 

The transverse momentum distribution for each cen- 
trality bin was corrected separately for the geometrical 
acceptance of the detector, the inefficiency of the tracking 
algorithm, secondary and incorrectly reconstructed par- 
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FIG. 3: Nuclear modification factor, Raa(pt), at mid-rapidity in bins of fractional cross-section at ^/s^ = 200 GeV (top row) 
and ^/s NN — 62.4 GeV (bottom row), for Cu+Cu (filled symbols) and Au+Au (open symbols). Systematic errors are shown 
with brackets (90% C.L.). The gray band represents the relative uncertainty on (N co u). 



tides, and the distortion due to binning and momentum 
resolution. The relative importance of these corrections 
and their estimated systematic uncertainties are similar 
to those reported in the previous analysis [f| . 

In Fig. we present the invariant yield of charged 
hadrons as a function of pr, obtained by averaging the 
yields of positive and negative hadrons. Data are shown 
for each centrality bin at both energies and are averaged 
over a pseudo-rapidity interval 0.2 < rj < 1.4. 

The centrality evolution of Raa{pt) at mid-rapidity 
for Cu+Cu collisions is shown in detail in Fig. [3J As 
a comparison, we also include the results from Au+Au 
collisions 0, |]| using the same centrality binning. As 
indicated, the most peripheral bin shown for 62.4 GeV 
Cu+Cu is 35-40% central, due to the limited efficiency 
of our vertex reconstruction for low multiplicity events. 

At both collision energies, we notice that, in the same 
bin of fractional cross-section, Raa is systematically 
higher in Cu+Cu compared to Au+Au. At y/s NN = 
200 GeV (top row), the Cu+Cu spectra exhibit a high-px 
suppression, relative to binary collision scaling, ranging 
from approximately 0.5 in the most central events to 
virtually no suppression in the most peripheral events. 

For both collision energies, the spectral shape of central 
Cu+Cu events appears to be very similar to the shape 
of peripheral Au+Au events at the same number of 
participants. This is illustrated in Fig. 0] where we 
present Raa versus the number of participating nucleons 
for both Cu+Cu and Au+Au collisions 0, § in bins of 
Pt- Since the previously published PHOBOS results did 
not achieve the statistics that now allow measurements of 
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FIG. 4: Nuclear modification factor, Raa, in bins of pr 
versus N par t at ^/s NN — 200 GeV, for Cu+Cu (filled symbols) 
and Au+Au (open symbols) 2, 4]- The gray band in the 
first frame represents the relative uncertainty on (N co u) and 
the solid lines show the effect of this uncertainty on the 
measured Raa- At high pr, bands are shown representing 
the predictions of a parton quenching model 47]. 

Cu+Cu spectra out to px — 7 GeV/c, we include results 
from the PHENIX collaboration 0] for comparison. 
The result illustrated in Fig. 0] is strikingly simple. 
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FIG. 5: Nuclear modification factor, Raa, in bins of pr versus 
Npart at y/s NN = 62.4 GeV, for Cu+Cu (filled symbols) and 
Au+Au (open symbols) 

Over the broad range of px that we measure, the bulk 
particle production seems to depend only on the size of 
the overlapping system, that is the Cu+Cu and Au+Au 
spectra are similar for the same number of participants 
(or binary collisions, see Fig.^). Although the measured 
centrality ranges in Au+Au and Cu+Cu collisions have 
less overlap at y/s NN = 62.4 GeV, this observation 
appears to hold at the lower energy as well (Fig.[SJ). 

At high-pr, a number of predictions have been made 
for the Cu+Cu system using models that successfully 
describe the centrality dependence of hadron yiel ds and 
back-to-back correlations in Au+Au 0, 0, UJ. One 
such model is the Parton Quenching Model (POM) [ft| . 
which utilizes BDMPS quenching weights [2(j and a 
realistic collision geometry to describe partonic energy 
loss. Although the centrality evolution of the Au+Au 
spectra are well fit by PQM, our Cu+Cu results suggest 



that this model slightly overestimates the suppression in 
the smaller system as shown in Fig. thus this predic- 
tion does not exhibit the observed N part scaling. Our 
conclusion is consistent with the prediction of a simple 
jet absorption model, whose only inputs are a Glauber- 
based collision geometry and a quadratic dependence of 
absorption on path length in an expanding medium |l9j ] . 
Using the absorption coefficient that describes the 200 
GeV Au+Au data, this model gives Raa values which 
depend only on N part , in agreement with our observation. 

Particle production at pr > 1 GeV/c in heavy- 
ion collisions is expected to be influenced by many 
effects. These include p^-broadening due to initial and 
final state multiple scattering (the 'Cronin effect'), the 
medium-induced energy loss of fast partons, and the 
effects of collective transverse velocity fields as well as 
parton recombination |2l| . Considering the significantly 
different geometries of Au+Au and Cu+Cu collisions 
with the same number of participant nuclcons, it is not 
obvious, a priori, that these effects should conspire to 
give similar spectra in both systems over such a large 
range in px- A full explanation of this observation, 
which appears to be a fundamental feature of heavy-ion 
collisions at these energies, may well present a challenge 
to theoretical models of heavy-ion collisions. 
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